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Surface Particle Motion of Rayleigh Waves in 
Prestressed Heterogeneous Orthotropic Elastic Half 

Space 

Inder Singh Gupta, Amit Kumar 


Abstract — The motion of the surface particles of Rayleigh 
waves in the prestressed heterogeneous elastic half space is 
discussed in detail. The analytical expressions for displacement 
components are derived and showed that these components are 
effected by initial stress and inhomogeneity factor. It is also 
proved that motion of the surface particles becomes retrograde 
elliptical in unstressed heterogeneous half space. For k=0, the 
ratio of the major and minor axes are derived and showed that 
the ratio becomes unity at the cut-off frequency. 

Index Terms — Rayleigh waves, prestressed heterogeneous 
medium, orthotropic medium, retrograde elliptic, retrograde 
circular.. 


I. Introduction 

The surface wave method is most widely used to find the 
mechanical properties of the medium, because it is 
non-destructive, has less inspection time, low cost and has 
wide range of applications in various fields like seismology, 
geophysics and materials science. 

The study of propagation of Rayleigh waves in 

heterogeneous half space have been studied among others 
Stonely(1934), Wilson(1942), Newlands(1950), 
Hook(1961), Dutta(1963), Karlson(1963), Singh(1965) and 
Sidhu(1970). A good amount of literature can be found in the 
standard books like pilant (1979) and A. Ben-Menahem and 
S .J. Singh (2000) etc. 

In most of above studies, the effect of inhomogenity factor in 
different forms are discussed for the study of propagation of 
Rayleigh waves. Though, earth is prestressed heterogeneous 
medium due to many factors like gravity field and variation in 
temperature etc. 

Biot (1965) was the first that has shown anisotropy was 
developed due to initial stress present in the medium. Using 
basic equation of Biot, some researchers like Sidhu and Singh 
(1983) hasshown the effect of prestress on the propagation of 
P, SV,SH seismic waves. 

Due to complexity of the problem, the study of transmission 
of Rayleigh waves in the prestressed heterogeneous half 
space was not enough. Das et al. (1992), abd-alla et al. (2009) 
and Kakar and Kakar (2013) have tried to drive the dispersive 
equationsof Rayleigh waves in prestressed heterogeneous 
media. Most of them used potential method to solve these 
problems. Norries (1983) also has pointed out that potential 

Inder Singh Gupta, Department of Mathematics, J.V.M.G.R.R. (P.G.) 
College, CharkhiDadri, Bhiwani 

Amit Kumar, Research Scholar, Department of Mathematics, Singhania 
University, Rajasthan 


N. 


rS 


c k \t$en 


Research Publication 


method is not suitable for prestressed media 

In the present paper authors want to discuss the motion of 
surface particles in the prestressed heterogeneous elastic half 
space using matrix method. The expressions for the 
components of displacement are derived. For that purpose the 
basic equations are taken from their previous paper (Gupta 
and Kumar, proceeding ISTAM, 2014). Here it is assumed 
that material properties and initial stress components are 
varying as 
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II. BASIC EQUATION 

Consider a semi-infinite, perfectly elastic prestressed, 
heterogeneous medium .The materials is either 
isotropic in finite strain or anisotropic with orthotropic 
symmetry .The principle directions of prestressed are 
chosen to coincide with the direction of elastic 
symmetry and the co-ordinate axes . The general 
equations of motion for prestressed solid in the absence 
of external forces are given by Gupta and Kumar 
(proceedingISTAM, 2014). 
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where B^ ± ,3^2 ^12 ^ 23 >Q± and 5-^,5^ are elastic 

coefficients and initial stresses in homogeneous 
orthotropic prestressed medium .(l°,|i°) ,p c are Lami’s 

constant and density of material at the free surface. 
Suppose the solutions of Eqs. (1-3) in the form of 
displacement components are as follow: 
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v — ( z) 8 ~ E x * -y”-* 3 1 - (4) 

w = 

where i = V - 1. Then eq. (4) define a plane harmonic 

wave propagating in the direction of the normal to the 
plane 4- yr\ = constant, (5) 

as shown in Fig. 1 
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Figure 1 
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On simplification of Eqs. (1-6) , we get 
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The set of Eq.(7) can be written in the following matrix form 

■y = at, (9) 

where T is column matrix with elements 

and A is 6x6 matrix(given by Gupta and Kumar) . Let the 

solution of Eq. (9)is of the form T=&* 2 C 

which is satisfied for all z if [A — sJ] C =0,(10) 
where I denoting six order unit matrix. 

Hh 

The set Eq.(10) has a non -vanishing solution vector C if and 
only if 

14-5/1 = 0 (11) 


The determinant of Eq.(ll)is a sixth order equation in s and 
gives the following six distinct values, 
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The general solution of Eq.(9) is Therefore 
f = £p =1 g *iz c 1 , (14) 

Writing the column vector I a no C , in the symbolic form as 

~ - _r ] , C t — [Cy], i,j - (1,2,3, ,6'i in Eq.(14) and 

equating the corresponding elements ,we get 

ui = Ti = Zf =1 ^ C a 

,»! = T 2 = If =1 s *i* C.j.uq = T 3 = Sf =1 Q a . 

(!5) 

We replace s by s s - bjic C iy C, in Eq. (10) ,we get 
jA - 5;/] [Cj] = 0 (16) 

Here we assume, that each of the sets of six C, : - for any given 

i n^dS#b e c ^ns'M eat^t means all six solution are L.D. ,we can 
express onefis^erms of L.C. of other five. For that, fixed i , 
and express them in such a way that none of the values of 

i=l,2,3 6, the coefficient determinate becomes zero. It will 

be possible, if we write Cg( j=2,3,....,6) as Qi^dsay). For 

each i, Eq. (16) gives a set of six homogeneous equations in 
six unknowns (j= 1,2, 3,...., 6)as 
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When the elastic half space is unstressed heterogeneuus 
infinite isotropic solidthen put P=R=0 in Eq. (8), we get 
l = = B° 2 = JL 0 + 2|£ 0 , = Q* = A, 

= m , y 3=(' r fli_y (21) 

□ ti \1-, h- 2 l! 


K 1 = 


Pa 


N. 


rS 


ext^en 

[Research Publication 


2 


www.ijntr.org 


International Journal of New Technology and Research (IJNTR) 
ISSN:2454-4116, Volume-2, Issue- 1, January 2016 Pages 01-05 


c = {X c + E M 2 - pV , 
D = (A 0 + 2ti 0 yr} 1 + 1 - p z p z , 

and 

E=p 0 a 2 - ftp 2 . 

After some manipulation, we obtain 
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Substituting for these in Eq.(15), we obtain 
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Appropriate solution to be used from Eq.(26) which satisfy 
the radiation as z -* go are f or K t - K 2 - A T 5 = 0 , where we 
must ensure that Re(.s Lf s 2f s z ) > 0. Therefore putting 
K 1 = K 2 = K 5 = 0 in Eq.(26),we get 
~ Sf= 3 ^L + K & e** z , 

^ = jEf =3 K- v - * , (27) 

vi- L =^Xf =3 

The boundary conditions are( Gupta and Kumar, proceeding 
ISTAM,2014) 
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Using Eq. (27) in Eq. (28), we get 

[«??+^/2)J 3 +((7|-/f 0 /2)tn 3 ]tf 3 +| (<3£+F°/2)s 4 + (Q£-S°/2) 

] 7f 4 + i +E " /2 )i e ,7 l £l =0 , 

f[(£’i — ^b)™s + 5i s i 

]*»+[(« i - S»)m 4 E 0i D **K*-4 (^(2^=0, 




(29) 


[(£12 E iijr E (^ja+-^gg) 1 ?"- 3 ' 2 j ?TC g s 3 ]^ 3 +[ 

(■^12 + ^iijf "+(^23 +^"g , g3 ; ! " — 32jjri 4 £ 4 ]A 4 =0 

If we pul <7 £ = QJ 1 ,5,", = 5g 3 and K & —0 then first two 

Eqs.Of Eqs. (29) becomes identical and from Eq. (29) , we get 

)m 1 ]K 1 +[Qls i + (Ql-S^ )m 4 ]K 4 = 0, 
ES&X 2 + (Sfg +j r “ 3 )fj 2 -Bi 2 m 2 s g ] Ag +[ 

(5£ : ES&X 2 +(3»+Sn>j a 
Elimination of and K 4 from Eq. (30) gives 

+ (Qj -5g° g )m g 

(B{ 2 E : + (Sfg E £33} I7 : - S-.VfTT.gSg 


(30) 


Q?s 4 + (Qi - 5°, }m 4 
(^12 E5»g)f ! + (S° +5«), J -S®m 4 s 4 


=0, 


(31) 


which is the frequency equation of Rayleigh waves 
propagating over the free surface of heterogeneous 
prestressed medium. 

On Simplification, the determinate value of Eq. (31) becomes 
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which give the displacements everywhere. For surface 
displacements we put z=0. FromEq.(30) 
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When the medium is heterogeneous unstressed (t. =0, 

R ' =0), then Eqs. (40), (42) and (43) reduce to 
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t is clear that the particles lie on the curve obtain from two 
3qs. (59) and (60). 

If we rotate the (u,v) axes to (u\. i/j through an angle of 0 
where cosO = - , sin0 =— then the relation between two 
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for s 3 and s 4 are real and positive iff 
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From Eq. (34) in this case 
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Where v is the direction of propagation of the Rayleigh 
waves and u is perpendicular to it. Then the curve 
represented by Eqs. (5 9) and (60) is equivalent to 


where LE is the value of D for E=0. 
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rju — <^v — 0 (59) 

and from first and third gives 


u w 

— ~ 4 2 — ~ = 1 

A*~ B 


(60) 


u — 0 


v 


A * 

rr 2A , 

a 'z2 



+ 


w 


B 



= 1 


(62) 


Using Eq.(61) and first two equations of Eq.(38) , we get 
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v = 


£ 


where A = a A 
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Eq. (62) shows that, the j(^Alcle paths are ellipses in vertical 
planes to the direction of Rayleigh wave propagation. 

Eqs. (63) and (64) gives 
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Hence Eq. (65), shows that the motion is elliptical. The ratio 
of the major and minor axes is given by Eqs. (57) and (58) is 
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Where the range of \| / is taken from the Eq. (55). Graphically 
we can see that the ratio of major and minor is increase with \| / 
as shown in figure (2). Again we have drawn a graph between 
- and — taking — along y-alis and — along x-axis for k=0. 
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IE. Stationary Rayleigh waves : 


For the surface displacements due to stationary Rayleigh 
waves we can take Eqs.(63) and (64) and change the sign of p 
and add, we get 

v = 2 A cos(x2 + yrj) cos pt , (67) 

w = 2B + yr / ) cos pt . (68) 


The particle motion is seen to be rectilinear simple harmonic 
and of the same frequency as that of the wave, if we put 


x4 + yr/ = cr 


( 4 *0 

x—+y— =crp l 
\ cr a) 


(69) 


where p ] is the distance from the origin in the direction of 


the wave propagation. Eqs.(67) and (68) can be written as 


v = 2 A cos(crp 1 ) cos pt , (70) 

w = 2Bsin( ap\) cos pt . (71) 
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From Eq. (6), the wavelength is K=— , then Eqs (67) and (68) 


give 

? 

, when 

K 3K 5K 

o 

ii 

> 

Pi a ’ A ’ A 

4 4 4 

A K 2K 3K 

ii 

o 

, when 

Pi-0’ 

2 2 2 


Hence, at interval of one fourth of wavelength , the 
disturbance is alternately vertical and horizontal. 


IV. Conclusion 

it is clear that the particles lie on the curve obtained by 
intersection of two surfaces represented by two Eqs. (59) and 
(60). From Eq. (65) we conclude that the particle paths are 
ellipses in vertical planes parallel to the direction of Rayleigh 
wave motion. The motion is retrograde elliptical if the signs 

of A and B are same for X =0. For different value of 


y/ lying in the range \ < y / < ^ + ^ ,we can draw the graph of 


B 

— in unstressed state of medium. When A =0, then iff =1 
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and — becomes unity at the cut-off frequency. Hence the 


particle motion is retrograde circular. Other than this value of 
the frequency , the particle paths are retrograde elliptical. It is 
clear from the figure (3) that the increasing value of frequency 
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upper most value, of the frequency, phase 


velocities approach Rayleigh wave velocities of the 
homogeneous case equal to 
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